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ABSTRACT 

We explore the halo structure of four gravitational lenses with well-observed, thin Einstein rings. We find 
that the gravitational potentials are well described by ellipsoidal density distributions in the sense that the 
best-fit nonellipsoidal models have parameters consistent with their ellipsoidal counterparts. We find upper 
limits on the standard parameters for the deviation from an ellipse of |af/arj| < 0.023, 0.019, 0.037, and 
0.035, and \a%/a \ < 0.034, 0.041, 0.051, and 0.064 for SDSS J0924+0219, HE 0435-1223, B 1938+666, 
and PG 1 1 15+080, respectively. We find that the lens galaxies are at the centers of their dark matter halos, and 
obtain upper limits for the offset of each center of mass from the center of light of | Ax| < 0."004, 0"005, 0"009, 
and 0."005, corresponding to 22, 29, 70, and 23 pc. These limits also exclude the possibility of any significant 
lopsidedness of the dark matter halos and set an upper limit of / sat < y/N% on the mass fraction of massive 
substructures inside the Einstein ring if they are divided over N satellites. We also explore the properties of 
galaxies as substructures in groups for the lens PG 1 1 15+080, finding evidence for dark matter halos associated 
with the galaxies but no evidence for a clear distinction between satellite and central galaxies. 
Subject headings: cosmology: observations — dark matter — gravitational lensing — quasars: individual 
(B 1938+666, HE 0435-1223, PG 1 1 15+080, SDSS J0924+0219) 



1. INTRODUCTION 

In cold dark matter (CDM) models, a typical virialized ob- 
ject consists of a luminous galaxy centered in a dark matter 
halo and surroun ded by satellites, some of which may con - 
tain no stars (e.g.. lBullock et al.ll2000l iKravtsov et all '2004). 
In this paper, we explore three questions about the structure 
of such halos using four gravitational lenses in which we can 
accurately measure the structure of the Einstein ring image of 
the quasar host galaxy. 

First, we explore the angular structure of the dark matter 
halos. The angular structure of luminous early-type galaxies 
is well-approximated by projected ellipsoidal di stributions, or 
ellipses, albeit with sm all systematic deviations (Bend eret al.l 
119891 iRest et aljf2001l) . Considerably less is known about 
the angular structur e of dark matter halos (e.g., lKatzlll99ll 
Allgo odet al.l l2005). particularly after their structure is mod- 
ified by the dissipative cooling and reheating of the baryons 
by star formation as the galaxy evolves. This issue is rele- 
vant to the halo structures less from an interest in the ellip- 
ticity of the halos than from the possibility that strong devia- 
tions from an ellipsoidal density distribution may be respon- 
sible for the flux ratio anomalies in gravitational lenses (e.g., 
| Mao & Schneider! fl998t IZhab & Pronkl l200ll IChibal 120021 
Metcalf & Zhao 2002) that otherwise provide the strongest 
evidence for the ex istence of dark matter sub st ructures in the 
halos of galaxies tDalal & Kochanekl 120021) . lEvans & Witt! 
(2003) demonstrated that lens potentials with arbitrary an- 
gular structure could reprod uce the flux a nomal i es without 
the su bstru ctures, but iKoc hanek & Dalall (120041) . lYoo et all 
(2005), and Congdon & Keeton (2005) demonstrated for sev- 
eral lenses that the deviations in the lens potentials from an el- 
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lipsoidal distribution are too small to remove the flux anoma- 
lies. Here we extend our study of the shape of the gravitational 
fields in galaxy halos to three more gravitational lens systems. 

Second, we explore whether the center of mass of the lu- 
minous galaxy is identical to the center of mass of the halo, 
since most popular halo models assume that luminous galax- 
ies are centered on t heir dark matter h alos (e.g.J Seliak 2000; 
iBerlind & W einberg 120021 ICoorav & Shethll2002l) . Although 
it is unlikely that the core of the stellar and dark matter dis- 
tributions can be offset in a galaxy, potentially there could be 
some lopsidedness or other effects that would make it invalid 
to assume that the center of light is the center of mass. Weak- 
ening the centering of the h alo on the lens galax y was also 
one of the factors permitting Evans & Witt (2003) to explain 
flux ratio anomalies with changes in the lens potential. Small 
numbers of massive, but dark substructures would also create 
an offset between the centers of mass and light. 

Finally, we explore the halo structure of a galaxy group. It is 
likely that central and satellite galaxies in a group show signif- 
icantly different mass-to-light ratios on large scales due to the 
mass differ ences betwee n halos associated with each galaxy 
(e.g., IZheng et al.ll2()0 5). This could explain, for example, 
the increasing evidence that e arly- type galaxies are hetero- 
geneous in their structure (e.g., Kochanek et al. 2005). Most 
lens galaxies are members of small galaxy groups, where the 
other galaxies in the group can perturb the lens through their 
tidal gravities. If we consider a lens centered on one of the 
satellite galaxies, then we might plausibly identify the central 
galaxy if one of the other galaxies in the group has a signifi- 
cantly higher mass-to-light ratio on large scales than the oth- 
ers. If, on the other hand, the structural differences between 
central and satellite galaxies are small compared to the intrin- 
sic scatter, then no such identification can be made. We can 
explore these possibilities using lens galaxies in small groups. 

We examine these three issues using four gravitational 
lenses with well-defined Einstein ring ima ges of the 
quasar host galaxy : SDSS J0924+02 1 9 (llnada et al J 1 2003 1 : 
' (X)5l). HE 0435-1223 dWisotzki et al.l 12000: 
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iMorgan et alJ l200l iKochanek et alJ l2005l). B 1938+666 
JKingetalJ Il99l~l1998t iTonrv & Kochanekl EOOOT) and 
PG 1 1 1 5+080 jW evmann et all 1198(3: ICourbin et all 
19971 iKeeton & Kochanekl Il997t llmpev et all Il998t 
Treu & Koonmans 2002). We use the shape of the Ein- 
stein ring formed by quasar host galaxies to accurately 
constrain the an gular structure of the underlying potential 
llKochaneket alJEOOTt IY00 et al] 12001 . In |2] we briefly 
present the Hubble Space Telescope {HST) observations 
of the four lenses followed by a short summary of the lens 
models in Our main results on the first two questions are 
presented in 50] We discuss the third question for the galaxy 
group surrounding PG 1115+080 in |5] We summarize our 
results in ^6] 

2. OBSERVATIONS 

We analyzed four gravitational lenses (SDSS J0924+0219, 
HE 0435-1223, B 1938+666, and PG 1115+080) with well- 
defined, thin Einstein ring images of the quasar host observed 
at //-band (NIC2/F160W) with the Near-Infrared Camera and 
Multi-Object Spectrograph on the HST. The data were re- 
duced and modeled usin g our standard procedures (iMcLeodl 
[T997t ILehar eUd]l2000h. and we extracted the Einstein ring 
curve as described in Kochaneket al. (2001). The image is 
fitted by a model consisting of the lens galaxy, the quasar im- 
ages, and a lensed host galaxy. The lens galaxy and quasar 
images are then subtracted, and we measure the radius r(9) 
of the Einstein ring as a function of position angle around the 
lens galaxy. The radius r{6) of the ring curve can then be used 
as a constraint on lens models under the assumption that the 
isophotes of the host galaxy are close to ellipsoidal. 

SDSS J0924+0219 consists of four lensed quasar images 
of a source at z, = 1 .524 produced by an elliptical g alaxy at 
zi = 0.393 JOfek et alJl2005l fEigenbro d"etai1l2005h and the 
lens is well- separat ed from nearby galaxies Jlnada et alJ2003t 
IKeeton et alJ2005l) . Eight dithered images of the Einstein ring 
were obtained on 2004 June 1 for a total integration time of 
4640 seconds (Morg an et al. in pre paration). We use the im- 
age positions from Keeton et aljl20 05). The lens is quite iso- 
lated and we use only an external shear to model its environ- 
ment. By isolated, we mean that no nearby object provides 
higher-order perturbations than an external shear. 

The lens system HE 0435-122 3 is composed of a four- 
image lensed quasar at z s = 1.689 ( Wisotzki et al. 2000) and 
an elliptical lens galaxy at zi = 0.4546 ( Morgan et al .1 120051 
lOfek et a"fl 120051) . A spiral-rich group of about 10 galaxies 
was found within 40" and the closest spiral galaxy G22 (SBb) 
is separated from the lens galaxy by 4 "46. We use the po- 
sitions of the lensed images and field galaxi es around the 
lens fr om the HST I 'ACS images described in Morgan et al. 
(120051) and the 2560 s HS77N ICMOS image describe d in 
Koch anek etall (120051) . As in iKochaneke t al. (200j), we 
model the environment with an external shear and the nearby 
galaxy G22. 

B 1938+666 co nsists of a normal early-typ e lens galaxy of 
redshift zi =0.88 (Tonrv & Kochanek 2000) that produces a 
two-image system and a four-image syst em of a radio source 
of unknown redshift (Kingetal. 1997). The host galaxy 
is seen as an alm ost perfectly circular ring around the lens 
jKing et al.lll998l) . We analyze the 2816s NIC2 image taken 
on 1997 October 7. Fo r the radio image positions, we use the 
5 GHz observations of King et al. ( 1997). We allow the reg- 
istration of radio and optical coordinates to be optimized as 
part of the fit. The lens is fairly isolated and we model the 



environment using only an external shear. 

PG 1 1 15+080 consists of four images of a z, = 1 .72 quasar 
lensed by a zi = .31 elliptical galaxy ( Wevmann et alJll980t 
Impe vet alJ I fl998l) . This system belongs to a small group of 
ga laxies producing non-negligible higher-order perturbations. 
In lYoo et alJ ll200 5). like most previous studies, we modeled 
the environment as a separate group halo. Here we explore 
mo dels with halos ce ntered on the galaxy positions provided 
bv llmpev et al] lll998l) . We use the image positions from 
Impe v etalJ (119981) . 

We use a concordance cosmological model with density pa- 
rameters fl„, = 0.3, ft\ = 0.7 and the Hubble constant h = 
Z/o/100 kms _I Mpc _I = 0.7 to calculate angular diameter dis- 
tances in a ACDM universe. 

3. LENS MODEL 

We use the scale-free potential with arbitrary angular struc- 
ture 

- cj)(r, 9) = F(9) = fl + a 2 cos 29 + b 2 sin 29 
r 

m=2 

00 

+ [ Aa m cosmfl + Ab m sin m9] , (1) 

m=3 

as our lens model. The model has a flat rotation 
curve, which also seems to be true of the typical lens 
galax y (e.g., iTreu & Koopmansl 120021 iRusin & Kochanekl 
2005). This lens model has been used extensively in 
other studies (IZhao & Pronkl I200H IE vans & Wittl 120011 
2001 IKochanek et all20pitlWucknitzl2002tlYoo et1illl2005t 
Congdon & Keeton 2005). The model parameters are ao, ai, 
b 2 , Aa m and Ab m with m > 3. The coefficients a\ m for the el- 
lipsoidal part of the surface density of the lens galaxy are de- 
termined from the quadrupole moment of the galaxy defined 
by the model parameters as 

3 

a\ = (fl2Cos20L + /72sin20£) , (2) 

fl 

and the major axis orientation, 9i = 0.5 tan" 1 (bijaj). Note that 
the m = 1 terms in this potential are degenerate with a shift in 
the source position and can be neglected. The standard param- 
eters for the deviation of the mass density of the lens galaxy 
from an ellipsoid are a„/flo = ( 1 -m 2 )Aa m /ao and b^/ap = 
(l-m 2 )Ab m /a() for m > 3 (e . g. . lBender"&M ollenhof fl 1 9871 
iBender et alJl!989T) . We consider either ellipsoidal models 
( Aa m = Ab m = 0) or nonellipsoidal models where Aa m ^ 
and Ab„, for m < 5. In either type of model, we expand 
the ellipsoid to order m = 5, beyond which higher-order coef- 
ficients and deviations are negligible for the relatively round 
lens galaxies (of axis ratio qi— 1) we consider here. 

We also include two sources of external perturbations. The 
first is an independent external shear d ue to objects distant 
from the lens or along the line-of-sight (Keeton et al. 1997). 
Thus, we add a term to the potential 

<Mr, 9) = i 7 r 2 cos 2(0 -0J = r 2 G 7 (0), (3) 

with two free parameters for the shear amplitude 7 and ori- 
entation 9 1 . The second source of perturbations are galax- 
ies near the main lens that can produce hig her-order pertur- 
bations beyond an external shear (e.g., Ke eton & Zabludo"fl 
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TABLE 1 
Best-Fit Models 



SDSS J0924+0219 



HE 0435-1223 



B 1938+666 



Parameter 


Ellipsoidal 


Nonellipsoidal 


Ellipsoidal 


Nonellipsoidal 


Ellipsoidal 


Nonellipsoidal 


Is 


0.86 ±0.06 


0.84±0.16 


0.47 ±0.02 


0.51 ±0.06 


0.56±0.10 


0.56±0.18 




OA A J- 1 A T. 


OS n 4- OO 1 
— Zo.U It zz. / 


— J.y ± Z.U 


— 0. 1 It J.Z 


1 4- O 3 
— 1 l.Z It V.J 


— 1 1 .Z It Z J. J 


7 


U.uJO 1 \J.U1 -J 


058 4- 023 


035 4- 003 






n riOfi 4- 007 


0^, ( deffs^ 
w t v xj ^'&' 3 / - ■ 


65.6±0.9 


62.0±3.5 


-60.2 ±5.1 


-68.3 ±4.1 


33.1 ± 14.0 


33.3 ± 15.6 


7g 






0.040 ±0.001 


0.042 ± 0.003 






r ? (arcsec) 


_ 


_ 


= 4.46 


= 4.46 


_ 


_ 


g (degs).. 


- 


- 


= -144.7 


= -144.7 


- 


- 


x/,(arcsec) 


-0.001 ±0.001 


-0.001 ±0.003 


-0.003 ± 0.002 


-0.001 ±0.003 


-0.007 ± 0.007 


-0.007 ± 0.007 


Vz,(arcsec) 


0.001 ±0.001 


-0.004 ± 0.003 


0.000 ±0.001 


0.001 ±0.003 


0.004 ±0.007 


0.004 ± 0.009 


HL 


0.87 ±0.01 


0.85 ±0.02 


0.62 ±0.01 


0.55±0.01 


0.89 ±0.02 


0.89 ±0.01 


e L (degs).. 


-68.1 ±7.3 


-67.2 ±9.8 


-14.3±9.1 


-10.5 ±9.4 


-54.9 ±8.5 


-54.5 ± 10.2 


C1Q 


0.88 ±0.01 


0.88 ±0.02 


1.21 ±0.01 


1.21 ±0.01 


0.46 ±0.01 


0.46 ±0.01 


10 2 a 2 


-1.4 ±0.2 


-1.6 ±0.4 


8.3 ±0.4 


11.1 ±0.5 


-0.3 ±0.3 


-0.3 ± 0.4 


10 2 b 2 


-1.4 ±0.4 


-1.7 ± 1.1 


^4.5 ±0.4 


-4.2 ±0.7 


-0.8 ±0.1 


-0.8 ±0.3 


10 3 Aa 3 .... 


= 


-0.9 ±2.5 


= 


2.1 ± 3.3 


= 


0.0 ±4.8 


10 3 A6 3 .... 


= 


-2.2 ±3.8 


= 


1.4±4.7 


= 


0.1 ±4.3 


10 3 Aa 4 .... 


= 


0.8 ±3.0 


= 


2.9±3.3 


= 


0.0±2.1 


10 3 A6 4 .... 


= 


1.2±3.3 


= 


0.1 ±6.6 


= 


0.0±1.7 


10 4 Aa 5 .... 


= 


-1.7 ± 11.9 


= 


0.1 ±8.4 


= 


0.1 ± 1.2 


10 4 Afo 5 .... 


= 


-1.1 ±15.6 


= 


0.1 ±5.3 


= 


0.0±1.3 


X 1 


54.9 


47.8 


20.9 


15.9 


4.1 


4.1 


x 2 


0.3 


0.1 


0.2 


0.2 


2.3 


2.3 


-Mens 
Xtol 


0.3 


1.8 


1.3 


0.4 


0.7 


0.7 


55.5 


49.6 


22.5 


16.4 


7.1 


7.1 


N<lof 


69 


63 


70 


64 


66 


60 


F-Test(%) 




93.4 




35.9 




70.4 



NOTE. — The best-fit ellipsoidal and nonellipsoidal models. Only astrometric constraints are used in the fits. Angles are standard position angles while the 
lens position and the coefficients of the lens potential are calculated in Cartesian coordinates in which A-direction is toward the West. The model parameters are 
the axis ratio q s and major-axis position angle 0„ of the source, the external shear (7, (9 7 ), the shear perturbation j g produced by a nearby galaxy at g ), 
the position of the lens (xfc, yl), the axis ratio qi and major-axis position angle 81 of the lens, the ellipsoidal lens potential {ao, ai, bi\, and the higher-order 
deviations Aa w and Ab m for m = 3, 4, 5 from an ellipsoid. 



120041 IKochanek et alJl2005l) . We approximate these galaxies 
by the potential, 



v (r,0) ~ r^aicosm(6-6 g ) + r 2 ^bf n cosm(6-6 g ) 



m=l 



m=\ 



^ env (0) + r 2 G env (0), 



(4) 



where the coefficients a% and b% m for m < 3 are determined by 
minimizing the difference between this model and the deflec- 
tion produced by a singular isothermal sphere (SIS), 0sis(f) = 
b g \r— r g \, where b g is the Einstein radius of the nearby galaxy 
and r g is its position vector from the lens center. The approx- 
imation is accurate to within 0"001 over an annulus of width 
1" encompassing the Einstein ring for the closest perturbing 
galaxy in our sample (galaxy G22 in HE 0435-1225, which 
is 4."5 from the lens galaxy). 

The total potential is 0(r,0) = r[F(0) + F env (0)] + 
f 2 [Genv(0) + G 7 (0)], which may include contributions 
from several nearby galaxies. For an ellipsoidal host galaxy 
with a monotonically decreasing surface brightness, the 
Einstein ring curve defined by the peak surface brightness of 
the Einstein ring as a function of angle going around the 
lens galaxy is simply 

hSt + u St 

^ )= — rs~ t — ' (5) 

where h = (F + F env )e r + (F' + Fl m )eg, the source plane tan- 
gent vector is t = M" 1 -e r = (l-2G en v-2G 7 )e r -(Ge nv +Gpee, 



M -1 is the inverse magnification tensor, the source center is 
«o, and S is the two-dimensional shape tensor of the source 
with axis ratio q s and position angle o f its major axis 0;, (see 
IKochanek et allioOltlYoo et all2005t) . 

For each lens we fit the image positions, the lens galaxy po- 
sition, and the Einstein ring curve using a simple x 2 statistic, 



x 2 (/» = £ 



|[Ho-n,-(p)]M,f ^ |r(0,)-r(0,;/7)| 2 



2 

"/.ring 



(6) 



where er,',i mg and airing are the uncertainties in the i-th im- 
age position and the ;-th Einstein ring radius, respectively. 
To allow rapid parameter searches, we use this source plane 
statistic for the \ 2 weighted by magnification M = |M| to ap- 
proximately correct for the difference between deviations on 
the source and image planes. Our fiducial ellipsoidal model 
has 1 1 free parameters, p; the source position is u = (u x , 
m v ), the axis ratio q s and major-axis position angle 9 S of the 
source, the ellipsoidal lens potential {ao, the posi- 

tion of the lens (x^, y£), and the external shear amplitude 
and direction (7, 7 ). In addition, there is one parameter 
7 g = b g /2r g for each nearby galaxy with observed position 
r g = (r g , Og) that may be added depending on the environment 
of the lens galaxy. Although flux ratios can provide additional 
constraints, they can be changed b y dust extinction (e.g., 
iFalco et aTJl 19991) . microlensing (e.g., Chiba et al. 20 05b, an d 
perturbations by substructure (e.g., Kochanek & Dalai 2004), 
so we will use only astrometric constraints. We used the 
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Levenberg-Marquardt method to minimize the \ 2 statistic of 
a given model and to estimate the parameter uncertainties, 
and we checked the convergence of our solutions by repeat- 
ing the minimization with a wide range of initial parameters 
and whether we have reached a global minimum by using the 
comp lementary downhill simplex method (e.g., iPress et aTl 
11992) . 

4. THE STRUCTURE OF LENS GALAXIES 

In this section, we analyze the lens systems to investigate 
whether the lens galaxies show any evidence for nonellip- 
soidal structures in their density distribution. As in Yoo et al. 
J2005I) . we first fit each lens as an ellipsoid centered on the 
center of brightness of the lens galaxy and then fit it allowing 
deviations from ellipsoidal structure or offset positions. The 
significance of the changes can be evaluated by using the F- 
test to estimate the significance of the improvements to the fit 
from adding the additional structural parameters and by con- 
sidering the scale of the deviations from ellipsoidal structure. 
The best-fit parameters and the F-test results for each model 
and each lens galaxy are presented in Tabled Each lens is 
modeled with the environmental components (shear, nearby 
galaxies) described i n j2 | We focus on higher order (m > 3) 
angular structure in 34. H and on th e ma tch between the center 
of mass and the center of light in 34.21 

4.1. Angular Structure 

Our basic result is that SDSS J0924+0219, HE 0435-1223 
and B 1938+666 appear to be ellips oids, a s we had previously 
found for PG 1 1 15+080 (lYoo et alJl2005l) . In all four cases, 
the improvement in the goodness of the fit from adding the 
m = 3 - 5 terms to the potential is statistically insignificant. 
The results of the fits are summarized in Tabled 

Figure ^ illustrates the differences in the critical lines 
and isodensity contours for the ellipsoidal and nonellipsoidal 
models of the three lenses. Figure|2]displays the values of the 
center of mass offsets and the higher order deviations bf n /a.Q 
(m = 3,4,5) from an ellipsoid. In Tabled we quantify the dif- 
ference of the two best-fit ellipsoidal and nonellipsoidal mod- 
els by providing the F-test probability for the nonellipsoidal 
model being consistent with the ellipsoidal model. In none of 
the lenses do the additional structural parameters of the nonel- 
lipsoidal model improve the fits. 

Although we generally ignore the constraints on the flux 
ratios in this paper, we note that the best-fit ellipsoidal mod- 
els reproduce the observed flux ratios except for the images 
with flux anomalies that are probably due to substructure or 
stellar microlensing. For example, the flux anomaly A1/A2 of 
PG 1115+0 80 turns out to be due to the stellar microlensing 
l lChiba et al.ll2005l) . The image D of SDSS J0924+0219 and 
the images A and C of HE 0435-1223 also show deviations of 
predicted compared to observed fluxes while the global lens 
models in both cases provide a reasonable fit to the image 
positions and the flux ratios of the rest of the images. Since 
the magnification is more sensitive to the higher order struc- 
ture in the potential, the best-fit nonellipsoidal models using 
only astrometric constraints fit the observed flux ratios of the 
anomalous images poorly. However, even if we impose the 
flux ratio constraints including the anomalies, the best-fit pa- 
rameters change little from those for the ellipsoidal model, 
and it is simply impossible to fit the flux ratios while simul- 
taneously obtaining a good fit to each Einstein ring. We con- 
sider models based only on the astrometric constraints in the 
rest of this paper. 



Ellipsoidal models fit the observed image positions and the 
Einstein ring of SDSS J0924+0219 well with xlot = °- 80 - 
While a slightly better fit with Xa f = 0-79 is found using the 
nonellipsoidal model, the F-test probability that the additional 
variables were unnecessary is 93.4%, so the astrometric con- 
straints from the Einstein ring strongly favor the ellipsoidal 
model for the lens galaxy of S PSS J0924+0219 over n onel- 
lipsoidal models. Contrary to lEigenbrod et alJ J2005I) . we 
have no difficulty finding a model that simultaneously fits the 
quasar images and t he Einstein r ing. T he difference probably 
arises because feigenbro d et alJ (12005) compared model pa- 
rameters from two separate fits to each constraint rather than 
fitting both constraints simultaneously. 

HE 0435-1223 is also well fitted by the ellipsoidal model 
with Xdof = 0.32 (see Fig.Q. In fact, we overfit the data and/or 
the uncertainties in the Einstein ring curve are overestimated 
by 76%. The ellipsoidal model with a nearby SIS and an ex- 
ternal shear provides a fit consistent wit h the recent me asure- 
ments of the image time delays ( Kochanek et al. 2005J), and 
the parameters of the best-fit nonellipsoidal model are con- 
sistent with those of the ellipsoidal model. In addition to the 
F-test indicating that the additional variables of this nonel- 
lipsoidal model are unnecessary, the best-fit nonellipsoidal 
model requires a high external shear of 7 = 0.065 compared 
to 0.035 for the ellipsoidal model. Considering that we have 
taken out the contribution from the nearby galaxy G22, it is 
unphysical to have external perturbations of ~7%. The high 
amplitude of the external shear in the nonellipsoidal model 
compensates for the higher order deviations of the lens poten- 
tial from an ellipsoid. 

The lens system B 1938+666 exhibits doubly lensed (open 
squares) and quadruply (open triangles) lensed radio images 
whose sources ar e two radio lobes rather than a central core 
iKing et al.ll997l) . Note that the double radio images are off- 
set from the Einstein ring. For this lens, we therefore in- 
clude two additional source positions for the radio images off- 
set from the center of the host galaxy producing the Einstein 
ring. The nearly circular Einstein ring rules out any signifi- 
cant higher order deviations in nonellipsoidal models, and the 
best-fit nonellipsoidal model is consistent with the ellipsoidal 
model (see Fig.|2j. 

4.2. Center of Dark Matter Halos 

In the previous subsection we constrained the position of 
the lens galaxy to agree with the position of the luminous 
galaxy seen in the HST images. This is the standard as- 
sumption of almost all lens models. In principle, the cen- 
ter of mass could differ from the center of light if the halo 
were genuinely offset from the stars, if the halo were sig- 
nificantly lopsided, or if there were any sufficiently massive 
dark substructures inside the Einstein ring. To test these 
possibilities, we drop the constraints on the lens position 
and repeat the optimization of the ellipsoidal models to es- 
timate the offset (xi, y£) of the center of mass from the lens 
galaxy. Physically it seems unlikely that the galaxy light can 
be significantly offset from the dark matter halo on galaxy 
scales, and this is borne out by these results. We find that 
the center of mass can be offset from the center of light by 
only A cm = 0."001 ± 0"003, 0"003 ± 0"004, 0"008 ± 0"009, 
and 0"001 ±0"002 for SDSS J0924+0219, HE 0435-1223, 
B 1938+666, and PG 1115+080, respectively. These corre- 
spond to physical offsets of A cm = 5.4 ± 16.3, 17.4±23.1, 
61.9±69.6, and 4.6±9.1 pc. The nonellipsoidal models only 
increase the uncertainties in the offsets by roughly a factor 
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position predicted by best-fit models, while the filled ellipse at the center shows the predicted position, ellipticity, and position angle of the source. The caustic, 
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source for these images. Note that the other source position is nearly identical to the center of the host galaxy (| A«| ~ 0"023). 
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FIG. 2. — Center of mass offsets (top left) and the higher-order deviation coefficients of the best-fit nonellipsoidal models. 



of 1.5 and provide upper limits consistent with the values of 
the offsets. For a typical dark matter fraction ins ide the Ein- 
stein ri ng of e ~ 30% ( Koon mans & Treull2003t iRusin et all 
2003b), the offset in the center of mass of the dark matter halo 
is related to the offset in the center of mass by Ad m ~ A cm /e, 
so the limits on the offset of the dark matter halo are some- 
what weaker. Note that this co-location of the stars and the 
halo may not hold in clusters. For example, models of the 
cluster lens SDSS J1004+41 12require an offset of > 5.7 kpc 
between the cD g alaxy and the center of the cluster potential 
(Oguri et al. 2004), and there is some evidence that cD galax- 
ies are not at rest with resp ect to the rest frame defined by the 
other cluster galaxies (e.g. . Ivan den Bosch et alJ2005t) . 

The offsets also test the lopsidedness of dark matter ha- 
los and the existence of massive dark substructures. The 
luminous galaxies show no evidence for lopsidedness, be- 
cause they are well modeled by ellipsoidal de Vaucouleur 
profiles. Since dark matter halos are also supported by ran- 
dom particle motions, they should also have difficulty sup- 
porting a lopsided density structure, and this is borne out by 
the tight limits on Ad m . Finally, the lack of an offset also 
sets a limit on the existence of massive dark substructures. 
Given a typical subhalo mass function, most massive sub- 
structures dominate the total mass contained in s ubstructure 
JKlvnin et alll999HMoore et al.ll9?9t rTavlor & Babu ll2004 
iGao et alJl2004t IZentner et alJl2005l) . If we add a single ob- 



ject whose mass is fraction / sat of the mass inside the Ein- 
stein ring at radius r sat for the lens galaxy, it would shift 
the center of mass by A cm = /sat^sat- Thus, we have a limit 
that / S at < A cm / (r 2 at ) 1//2 where (r 2 ^) 1 / 2 is the rms distance 
of the satellites inside the Einstein ring from the lens cen- 
ter. Crudely, {'f at ) 1 ' 2 — b, the Einstein radius, so / sat < 1%. 
Adding additional satellites only weakly affects the limits, in- 
creasing the limit on / sat by \/~N for N equal mass satellites. 
This strong sensitivity to massive substructures is presumably 
one of the reasons that it is almost impossible to model lenses 
like MG J 04 14+05 34 wit h visible satellites using only a cen- 
tral potential fe.g.. lRos et all200CtlTrotter et alJ20 00). 

5. GROUP VS. GALAXY HALOS IN PG 1 1 15+080 

In the halo model that describes the relation of galaxies and 
dark matter, a luminous central galaxy is located at the center 
of dark matter halo and the satelli te galaxies are distr ibuted 
around the central g alaxy (e.g. , Berli nd & Weinberg! I2002t 
iKravtsov et alJ 120041) . In 34.21 we found that visible lens 
galaxies must be centered on their halos, and that they gen- 
erally lack massive substructures. If, however, we consider a 
lens consisting of a galaxy in a group, like PG 1115+080, 
then we can explore the relationship between central and 
satellite galaxies. This is potentially import ant because the 
time delay measurements for PG 1115+080 JSchechter et alJ 
H997I) require a main lens galaxy with almost no dark mat- 
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ter at the concor dance value of th e Hubble constant Ho = 
72 kms" 1 Mpc" 1 (iKochanekl 2003), in marked contrast to 
the typical lens galaxies, which seem to have a significant 
dark h alo (e.g.. iTreu & Koopmansll2002l: iRusin & Kochanekl 
2005). This case suggests that the structures of early-type 
galaxies are heterogeneous. One possible interpretation of 
this heterogeneity is that the primary lens of PG 1115+080 
is a satellite lens galaxy of the group that has been stripped 
of most of its dark matter. In the context of the halo model, 
this implies that one of the other visible galaxies around 
PG 1115+080 must be the central galaxy of the group, and 
it should show a significantly higher mass-to-light ratio than 
the other galaxies, all of which are satellites. In this section, 
we attempt to interpret PG 1115+080 within this theoretical 
picture. 

5.1. The Environment ofPG 1115+080 

Figure [5] shows the environment of PG 1115+080. Ear- 
lier models with an independent external shear or an SIS 
external perturber consistently indicate that the lens has a 
~10% shear perturbation associated with its parent group, 
with weak evidence that the group halo should be lo- 
cated nea r the luminosity weighted position of the nearby 
galaxies tS chechter et al j [1^27i iKeeton & Kochanek ll^^t 
I Courbin et al.lll997t llnroev et alii 19981 IZhao & Pronkll2001t 
Koch anek & Dalalll2004: lYoo et al.ll2005l) . The galaxies Gl, 
G2 and G3 are comparable in luminosity to the main lens 
galaxy, and the less luminous galaxy G14 is the closest. To 
estimate the perturbations produced by each galaxy, we use 
the /-band luminosities L of the nearby galaxies in Table |2] 
( Impe v~etalJl99 8) and the Faber- Jackson relation to estimate 
their stellar velocity dispersions cr+ oc L 1 / 4 and the Einstein 
radii b oc er 2 (L oc a\ oc b 2 ). Scaling all the galaxies based on 
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FIG. 3. — The environment of PG 1115+080. The nearby galaxies are rep- 
resented by the labeled circles, where the sizes of the circles are proportional 
to the perturbation produced by each galaxy. Solid circles are for the exter- 
nal shear, and filled circles are for the higher-order perturbations. The dotted 
concentric circles centered on the main lens galaxy are reference circles with 
radii of 10'.' 0, 20"0, and 30"0. 



TABLE 2 

Galaxy Environment of PG 1 1 15+080 



Object 


/(mag) 


r g (arcsec) 


# g (degrees) 


Is 


Lens 


18.55 ±0.50 


= 






Gl 


17.85±0.01 


23.285 


241.80 


0.030 


G2 


18.73 ±0.04 


11.957 


266.05 


0.039 


G3 


19.44 ±0.02 


18.543 


228.97 


0.018 


G8 


21.52 ±0.08 


30.624 


246.52 


0.004 


G9 


19.44 ±0.02 


14.352 


252.57 


0.007 


Gil 


22.98 ±0.34 


11.881 


252.62 


0.006 


G13 


23.70 ±0.10 


8.214 


113.16 


0.006 


G14 


23.30±0.10 


7.313 


154.82 


0.008 


G15 


23.80±0.10 


12.198 


167.09 


0.004 



NOTE. — Galaxy positions are relative to the lens. We use the /-band 
luminosities and the Faber-Jackson relation to estimate the relative am- 
plitudes of the external shear 7 g produced by each galaxy. The nonlinear 
perturbations are smaller than 7 ? by r g . The positions and luminosities 
are from Impev et al. 1998). 



the luminosity and Einstein radius of the main lens galaxy, 
we then compute the relative amplitudes of the external shear 
oc b/r g and the nonlinear perturbation oc b/r 2 produced by 
each galaxy given its distance r g from the lens (see Table|2ji. 
As illustrated in Figure|5J the luminous galaxies Gl, G2, and 
G3 dominate the perturbations while the nearby, less lumi- 
nous galaxies G13 and G14 make a non-negligible contribu- 
tion. 

We first test whether putting mass only at the positions 
of the visible galaxies can adequately fit the data, initially 
with no extra external shear. The galaxy positions are mea- 
sured to such high accuracy that we fix the observed position 
r g = {r g , 8 g ) for each galaxy, while we vary their Einstein radii 
bj through the free parameter j g We started with a model 
using just the three bright galaxies Gl, G2, and G3. Figure|4] 
shows the best-fit models. The three-galaxy environment fits 
the data quite poorly compared to standard models with a 
group halo centered near the center of light of the galaxies 
(A% 2 ~ 28). Only if we allow a large, independent external 
shear (7 > 0.07) can we obtain a good fit. However, when we 
include the next leading perturber, G14, the goodness of fit is 
somewhat better (Ax 2 — -5) than for the fit with a halo unas- 
sociated with the visible galaxies. Note however that this new 
model has one more degree of freedom compared to the stan- 
dard model, so the improvement in the fit is physically inter- 
esting but not statistically significant. Including the fifth most 
important galaxy, G13, hardly changes the goodness-of-fit, 
and the Einstein radii of the galaxies change little (although 
there are mass degeneracies between Gl and G2/G3). If we 
add an independent external shear to the four or five galaxy 
models, the estimated shear amplitude is small (7 ~ 0.01), 
and the goodness-of-fit and the Einstein radii of the galaxies 
remain unchanged. The addition of any more nearby galax- 
ies simply contributes to the external shear 7 because they 
produce negligible higher-order perturbations. Therefore, we 
adopt the four-galaxy model of Gl, G2, G3, and G14 with an 
external shear of amplitude 7 < 0.01 as our fiducial model for 
the environment in PG 1115+080. Next we address whether 
the masses of these galaxies show any correlations that test 
the structures of their halos. 

5.2. Halo Structure of the Galaxy Group ofPG 1115+080 

Our parameter for each galaxy, j g oc M(< r g )/r 2 , is propor- 
tional to the projected mass M(< r g ) of each galaxy out to the 
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distance of the galaxies r g from the lens. Based on the halo 
model, our working hypothesis is that one of these four galax- 
ies is the central galaxy of the halo and should have different 
halo properties than the other satellite galaxies. 

For any fixed assumption about the halo structure and as- 
suming the galaxies have similar mass-to-light ratios, we 
would expect to see a correlation between these mass esti- 
mates bi and the galaxy luminosities L,-. If there is a cen- 
tral galaxy with a very different halo mass than the satellite 
galaxies, then it should be identifiable as the galaxy least in 
agreement with the correlation of the other galaxies. We first 
consider two limiting cases in which the halos are either all 
truncated on a scale less than the distance to the lens or where 
all halos conspire to have flat rotation curves out to the dis- 

1 11 

tance to the lens. The Einstein radius is bf l oc L/ in both 
cases while the shear amplitude produced by each galaxy is 
7i = b 2 /r 2 gj in the former case and 7; = b,-/2r w - in the latter 
case. We include the correlation by adding a term to the fit 
statistic of the form 

^ \\nbi-\nbf\L,)\ 2 
A X 2 M . L = 2(N g . di - 1) lna + £ i ' ;J , (7) 

1=1 

where a is the logarithmic scatter in the correlation and N„d = 
4. 

Figure |5] compares the two mass models with and without 
an external shear 7. When we include the shear, we use a 
prior constraint of 7 = 0.00 ± 0.01 (based on the previous four 
galaxy model). The model assuming extended halos fits the 
data slightly better (A% 2 — 2) than the point-mass model re- 
gardless of the assumed scatter or the addition of the external 
shear. The logarithmic scatter of the correlation is a ~ 0.5 
for the best-fit model, and will give ^65% fractional dif- 
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FIG. 4. — Best-fit models for the environment in PG 1115+080. The 
goodness-of-fit relative to the best-fit model with a group halo unassociated 
with the visible galaxies is shown with the Einstein ring radii of galaxies in 
each model. The environment is approximated as three (open squares), four 
(open circles), and five nearby galaxies (filled circles) from top to bottom. 
Note that the models with four and five nearby galaxies are degenerate. 
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FIG. 5. — Mass models for the galaxy environment. We assume a Faber- 
Jackson relation L oc b with a natural logarithmic scatter <r. We show the 
goodness-of-fit as a function of the scatter for each mass model relative to 
the best-fit model of the four nearby galaxies with arbitrary mass. Solid lines 
show the results for galaxies with extended halos, and dashed lines show the 
results for point-mass halos. The light curves allow no additional external 
shear while the thick lines allow a small external shear. Since the mass scales 
of the galaxies are not in perfect agreement with the correlation in the mass 
models without an external shear, the difference in \ 2 becomes substantial 
as a tight correlation is imposed (cr <C 1). The mass models with an external 
shear can perfectly agree with the correlation, and the A\~ is controlled by 
the assumed prior on the scatter. Two vertical lines show the logarithmic 
scatters found in the SDSS for nearby elliptical galaxies 1 Bernardi et al. 2003) 
and for ensembles of lens galaxies 1 Rusin et al. 2003a). 



ferences. Note that there is little evidence of a useful cor- 
relation, but the scatter is not significantly different from the 
Faber-Jackson relation between the luminosity L and the cen- 
tral velocity dispersion cr* observed for early-type galax ies 
in the SDSS (a ~ 0.34 in our units, iBernardi et aill2003l) or 
that measured from ensembles of lens galaxies (cr ~ 0.24, 
Rusi neTaIl l2003a). When an external shear is allowed, it 
turns out that there is degeneracy between the external shear 
and the mass models for the four galaxies such that all the 
galaxies can perfectly agree with the assumed correlation, and 
the results are controlled by our prior on the scatter. The ex- 
tended halo model requires 7 ~ 2% to fit the data with all 
galaxies on the correlation, while the point-mass model needs 
7~4%. 

As our next experiment, we allow one of the four galaxies to 
have an arbitrary mass while imposing the correlation on the 
rest of the three galaxies (iVg a i = 3). Unfortunately, no single 
galaxy stands out, and the best-fit models without an external 
shear again require large intrinsic scatters cr ~ 0.5. When an 
external shear is allowed in the model, the external shear can 
always compensate for the shear produced by the galaxies due 
to the degeneracy between the external shear and the mass 
models, and no models are significantly better than the models 
with the four galaxies in the correlation. 

6. SUMMARY 
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In this paper we explored the structure of halos using four 
gravitational lenses with well-defined Einstein ring images of 
the quasar host galaxies. Adopting a scale-free potential with 
arbitrary angular structure, we considered both ellipsoidal and 
nonellipsoidal gravitational potentials, searching for any devi- 
ation from an ellipsoidal potential. Based on the astrometric 
constraints, none of the additional free parameters represents 
an improvement over the best-fit ellipsoidal models, and the 
best-fit nonellipsoidal models are consistent with their ellip- 
soidal counterparts. Furthermore, imposing the observed flux 
ratios as constraints does not alter the results. Anomalous flux 
ratios can not be eliminated by changing the central potential. 
The angular structures of all four lens systems are consistent 
with ellipsoidal models. 

We also investigated the offset of the center of mass of the 
lenses from the measured center of light by dropping the con- 
straints that the lens is centered on the visible galaxy. We 
find that the offsets are required to be A cm < 10 pc. This 
roughly implies that the dark matter halos can be offset by 
at most Ad m ~ A cm /e < 30 pc, where e ~ 30% is a typical 
fraction of dark matter inside the Einstein ring. This may be 
markedly different from clusters in which the cD galaxy can 
be significantly offset from the center of the cluster poten- 
tial. Dark substructures of mass fraction / sat could also shift 
the center of mass if there are few of them, and this allows 
us to set a limit of / sat < \f~N% for N equal mass substruc- 
tures inside the Einstein ring. This limit is not tight enough 
to represent a conflict with either theoretica l (iTavlor & Babul 
120041: iGao et al J 12004 iZentner et all 12005ft or observational 
(Dalai & Kochanek 2002) estimates of the substructure frac- 
tion. 



Finally, we explored the environment of PG 1115+080 to 
study the halo structure of its parent group under the assump- 
tion that the individual dark matter halos are centered on the 
nearby luminous galaxies. We focused on the four galaxies 
Gl, G2, G3, and G14 producing the largest perturbations, 
since the other nearby galaxies produce negligible higher- 
order perturbations and have effects degenerate with an in- 
dependent external shear. Placing the mass solely on these 
galaxies can fit the data, but we find only weak evidence for 
any systematic correlation between the mass and the lumi- 
nosity. Models assuming the galaxies have extended halos 
on the scale of their separations are slightly better than those 
where they do not, but the scatter in the correlation between 
the Einstein radius and /-band luminosity must be ~ 65%. 
This is comparable to that implied by Faber-Jackson relation 
observed for nearby galaxies or for lens galaxies. We also 
found no evidence for a marked difference between one of 
the galaxies (the "group halo central galaxy") and the others 
(satellite galaxies). The rapid convergence of the perturba- 
tions from nearby galaxies to a single external shear probably 
means that lenses like PG 1115+080 provide an insufficient 
number of constraints to test the hypothesis of the halo model 
in detail. 
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